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Abstract 

We define auto- and cross-correlation functions capable to capture dynamical characteristics induced by 
local phase space structures in a general dynamical system. These correlation functions are calculated in 
the Standard Map for a range of values of the non-linearity parameter k. Using a model of non-interacting 
particles, each evolving according to the same Standard Map dynamics and located initially at specific 
phase space regions, we show that for 0.6 < k < 1.2 long-range cross-correlations emerge. They occur 
as an ensemble property of particle trajectories by an appropriate choice of the phase space cells used in 
the statistical averaging. In this region of A:-values the single particle phase space is either dominated by 
local chaos (k < kc with kc ~ 0.97) or it is characterized by the transition from local to global chaos 
{kc < k < 1.2). Introducing suitable symbolic dynamics we demonstrate that the emergence of long-range 
cross-correlations can be attributed to the existence of an effective intermittent dynamics in specific regions 
of the phase space. Our findings further support the recently found relation of intermittent dynamics with 
the occurrence of cross-correlations (F.K. Diakonos, A.K. Karlis, and P. Schmelcher, Europhys. Lett. 105 , 
26004 (2014)) in simple one-dimensional intermittent maps, suggesting its validity for two-dimensional 
Hamiltonian systems. 
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I. INTRODUCTION 


Dynamical systems with mixed phase spaee (PS) [HI are eharaeterized by the eoexistenee of 
time seales whieh may differ by several orders of magnitude. The fast propagation in the ehaotie 
sea is usually interrupted by long-time stiekiness on loeal PS struetures sueh as dynamieal traps 
of hierarehioally arranged islands, net-traps or dynamieal traps of stoehastie layers [|2l. Usually 
the dynamies in these regions is almost regular in eontrast to the stoehastie profile of the eomplete 
trajeetory. Furthermore the PS geometry in the regions where regular and ehaotie trajeetories 
approaeh eaeh other arbitrarily elose is typieally fraetal and this leads to the presenee of unusual 
statistieal properties in the ensemble of the ehaotie trajeetories (Sl. This is mainly due to the 
faet that the dynamies is in this ease pseudo-ergodie [|4|| making impraetieal the use of trajeetory 
based eommon tools, like for example the Lyapunov exponents or the eorrelation funetions. In 
the eontext of Lyapunov exponents the main interest for a non-hyperbolie Hamiltonian system is 
to eapture the large variation of the loeal instability due to the eoexistenee of ehaotie and laminar 
phases along a referenee ehaotie trajeetory [|51. A suitable tool for this purpose is provided by 
the finite-time Lyapunov exponents [|6l whieh turn out to be mueh more sensitive to loeal PS 
struetures than the global Lyapunov exponents and allow for a better understanding of the mixed 
PS dynamies. This is elearly demonstrated in [|3 using as example Chirikov’s Standard Map 
(SM) [[81 whieh is the prototype of a Hamiltonian system with mixed phase spaee. It provides 
an approximate deseription of several physieal systems like the kieked rotor Q, the relativistie 
eyelotron [fTOl and the equilibrium eonfigurations of the Frenkel-Kontorova model [HTI . Being 
simple to simulate, two-dimensional and diserete in time, the SM has been extensively studied 
in the last deeades with emphasis on diverse dynamieal aspeets sueh as anomalous diffusion and 
stiekiness IfT^ . aeeelerating modes [[T3l and generating phase spaee partitions [fldl . whieh are 
elosely related to the presenee of loeal PS struetures and their impaet on mixed PS dynamies. 

For the ease of eorrelation funetions the implieations of pseudo-ergodieity are more intrieate, 
sinee these quantities eontain an additional time seale (the delay time) whieh implies a lower eut- 
off on the trajeetory lengths used in the eorresponding time-averaging. In faet for the eorrelation 
funetion to be representative for the PS dynamies it is required that the assoeiated trajeetory length 
is order(s) of magnitude larger than the maximum delay time. Thus the loeality of the PS struetures 
and the eonditions of long time propagation are eompelling faetors making the use of eorrelation 
funetions for exploring loeal PS dynamies a diffieult task. On the other hand, in a elass of Hamil- 
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tonian systems with PS dominated by stratifieation due to the presenee of invariant KAM spanning 
eurves, different PS regions may be dynamieally diseonneeted. Thus, ergodieity may be partially 
restored within an isolated PS region and the definition of eorrelation funetions within sueh a re¬ 
gion is plausible. Usually, in this ease, even after the destruetion of the invariant spanning eurves, 
remnants of the stratified strueture survive in the form of a dense set of stability islands [|T5l . Con¬ 
sequently, the stiekiness on this net of dynamieal traps induees large differenees between the time 
seale eharaeteristie for an ergodie eovering of a single PS zone through a ehaotie trajeetory and the 
time seale needed to eross the destroyed spanning eurve and enter in the aeeessible neighboring 
PS zone. One eould naturally think that even in this ease, despite the presenee of global ehaos, the 
definition of appropriate eorrelation funetions sensitive to the loeal PS struetures may be possible. 

The aim of the present paper is twofold. Firstly we introduee a elass of eorrelation funetions 
possessing the ability to explore the dynamies in restrieted PS regions for an arbitrary Hamiltonian 
system. Seeondly, using these eorrelation funetions we analyse the dynamies of the SM elose to 
the transition from loeal to global ehaos. In partieular we foeus on eross-eorrelation funetions 
between trajeetories with different initial eonditions. Sueh a trajeetory ensemble eorresponds to 
a set of non-interaeting partieles eaeh evolving aeeording to the SM dynamies with the same k. 
We show that around the loeal to global ehaos transition point, long-range eross-eorrelations de¬ 
velop, and subsequently, inspeeting the trajeetories in the ensemble used for the eross-eorrelation 
funetion ealeulations, we reveal the presenee of effeetive intermitteney in a symbolie representa¬ 
tion of the assoeiated dynamies. The observed intermitteney is strong in the sense that the mean 
waiting time in the laminar region diverges. Thus, our analysis demonstrates that the previously 
reeorded emergenee of eross-eorrelations in a system of non-interaeting partieles following non- 
Hamiltonian 1-d intermittent dynamies [[T^ is present also in a 2-D Hamiltonian system like the 
SM elose to the transition from loeal to global ehaos. 

Our paper is organized as follows. In Seetion II we give the definitions of the eorrelation 
funetions and we introduee the dynamieal system we are eonsidering in the subsequent analysis. 
Seetion III provides the numerieal results of the eorrelation funetions in the eonsidered model for 
a dense set of values of the non-linearity parameter k around the eritieal value kc signaling the 
transition from loeal to global ehaos. In Seetion IV we analyse the obtained results by performing 
a detailed evaluation of the eharaeteristie trajeetories and introdueing an appropriate symbolie dy¬ 
namies to reveal the underlying effeetive intermitteney. Finally Seetion V provides our eoneluding 
remarks. 
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II. CORRELATIONS IN THE STANDARD MAP - MODEL AND BASIC OBSERVABLES 

The model we consider consists of M non-interacting particles each evolving according to the 
dynamics of the Standard Map ([|8ll): 


pi+i sin 

= »!•’ + Pi‘1. 

where is the control parameter of the non-linearity and i = 1,2, We will here consider 

exclusively the case = . = k which simplifies the description significantly, 

reducing the study of the M-particle system to the study of an ensemble of initial conditions for 
a single particle. In this case the particle index [i] can be omitted in the description. Then the 
variables of the system, and stand for the single particle angular position and momentum at 
time instant n respectively, both variables calculated with modulo 27r. We restrict our analysis to 
the range k G [0, 2]. The system undergoes a transition from local to global chaos at the critical 
value kc = 0.971635... IfTTl . At kc the last KAM curve with ratio 0 (golden ratio) is destroyed 
and therefore chaotic orbits are not trapped between KAM tori any more, but they evolve across 
their remnants covering in principle all the available PS except the remaining islands of stability. 
Dynamically the trace of the destroyed regular PS structure is revealed in the stickiness of the 
chaotic orbits around these stability islands. As mentioned in the introduction it is likely to assume 
that the sticky evolution of the chaotic orbits in the immediate neighborhood of the remaining 
stability islands may introduce some correlations even between chaotic orbits originating from 
different initial conditions. To explore this scenario we will calculate the auto-correlation (AC) 
and cross-correlation (CC) functions of the considered system for the aforementioned fc-values, 
focusing mainly on the fc-region where the transition from local to global chaos takes place. 

For this task one usually employs the standard definition [[Tl for the correlation functions: 


Crii(m) = lim 

> N^OO 


^ N—m—1 ^ N—m—1 N—m—1 

— Xi{n)xj{n + m) - . Xi{n) 


N — m 


n=0 


(N — mY ^ 

^ > n=0 


n=0 


2 ) 


where N is the orbit’s length, x represents the PS variables 9 oxp,m stands for the usual time-delay 
occurring in the correlation functions (CFs) and i, j label the different sets of initial conditions 
(9q,po) that define the orbits i and j. The diagonal entries i = j in Eq. Q determine the AC 
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function which expresses the statistieal similarity of a single trajeetory at any two time instanees 
differing by m. The CC funetion is obtained from the non-diagonal entries i ^ j m Eq. ([^ and 
expresses the respeetive statistieal similarity between two different orbits (i.e. two orbits with 
different initial eonditions) at any two time instanees differing by m. The limit iV —)■ cx) is usually 
replaeed by averaging over the entire PS requiring the system to be strongly ergodie, whieh is not 
for all the eonsidered values of k the ease for the SM. In faet, in the most interesting region of k 
G [0.6, 2.0], the PS is mixed. As a result the assoeiated averaging eontains qualitatively different 
dynamieal eomponents (ehaotie, regular trajeetories) and washes out properties (i.e. eorrelations) 
whieh eharaeterize eaeh dynamieal eomponent separately. Thus, it is obvious that global averaging 
is not representative of the loeal struetures of the phase spaee. 

Another property that disfavors the use of the global averaging is related to the symmetries of 
the SM. The PS of the Standard Map is point-symmetrie around (vr, tt) eausing the faet that for 
eaeh trajeetory, denoted as (6*,p), there exists a partner trajeetory {9,p)ps = (27r — 6, 2tt — p) 
(where the index ps is used to notiee the point-symmetrie partner of {9,p)). This implies that the 
eontribution to the CC funetion of a pair of orbits {a,b} and the eontribution of the pair {a, bps} 
with bps the point-symmetrie partner of b, will be exaetly opposite eaneeling out in the averaging. 
Similarly, the eontributions to the AC funetion of a and Ups will also eaneel out in the averaging 
having equal magnitude and opposite sign. Thus, as long as there eoexist orbits in the ensemble 
with their point-symmetrie eounterparts, the averaging will always lead to vanishing auto- and 
eross-eorrelations. 

In order to overeome both, the peeuliarities of the averaging proeedure originating from the 
faet that the ergodieity is not strongly satisfied, as well as the trivial behavior of the eorrelation 
funetions due to the PS anti-symmetry, we introduee here a new elass of eorrelation funetions, 
the loealized finite-time eorrelation funetions (LFTCFs). Their aim is to extraet information on 
the possible emergenee of eorrelations due to the loeal struetures in the phase spaee and they are 
defined as follows: 


LC}'^\m) = lim — 

Sd^oo Sd ^ 

N—m—1 


where 


N — m 




xf\n)x^p{n + m) — 


N—m—l 


TV—m—1 


n=0 


(N-mY ^ ^ ^ xf\n + m) 

^ ^ n=0 


n=0 
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The LFTCFs denoted by LC in Eq. Q are obtained by the usual definition of CFs (Eq. Q) 
by fixing the trajeetory length N in the ensemble to a finite value and averaging over different 
CEs, each calculated using trajectory pairs with initial conditions within a specific PS cell 
The index d is used to indicate the PS domain where the cell belongs to (see below). In a 
similar manner as in Eq. ([^ the case i = j corresponds to the localized finite-time auto-correlation 
function (EETACE) and i ^ j corresponds to the localized finite-time cross-correlation function 
(EETCCE). With Sd we notice the total number of pairs of initial conditions used for the calculation 
of the CEs. In Eq. ([^ should necessarily belong only to one of the two distinguishable domains 

covering the entire PS: 

• The chaotic domain (d = c) which contains all the PS areas where local or global chaos is 
present. Eor k > kc the chaotic domain spans over the whole PS, excluding the stability 
islands. 

• The regular domain (d = r) which contains all the areas of the PS that are either stability 
islands or spanning curves (periodic or quasi-periodic orbits). 

As an example in Pig. we show the PS of the SM for k = 0.95. The chaotic domain corresponds 
to the colored regions while the regular domain to the black ones. 

Since the definition of the EPTCPs requires the averaging over initial conditions in a specific 
PS-cell one can naturally ask about the dependence of the form of these CPs on the specific loca¬ 
tion or size of the cell within a domain as well as the length of the involved trajectories. This will 
be thoroughly discussed in section III. Obviously the EPTCPs allow us to analyze the correlations 
in each dynamical component (chaotic, regular) separately, taking also into account the impact of 
the local PS structure. However the cell choice requires some care since one has to ensure that it 
belongs entirely only to one of the two aforementioned domains. 


III. NUMERICAL SIMULATIONS 

Eor finite Sd Eq. Q provides us with an estimate of the EETCCE with initial conditions in the 
cell Of course the result will in general depend on both the finite time interval N and the 
position of the cell This holds also for the EPTACPs. Eor a regular domain {d = r) the 
dominant frequencies define a characteristic time scale and N can be chosen to be a multiple of 
the period corresponding to the largest frequency peak in the related power spectrum. In addition 
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the regular motion does not generate diffusion and therefore the evolution of the initial eell 
remains localized in PS. Thus the validity of the terms finite time” and ’local” used in Eq. Q is 
straightforward for an ensemble of trajectories with initial conditions in a regular cell. 

In contrary, the validity of these terms is less clear when the cell of initial conditions lies within 
a chaotic domain {d = c). To illuminate this issue we explore in more detail the PS structure. 
For 0.55 < k < 0.85 local chaos appears in the vicinity of the large regular island centered 
around (tt, 0). For k E [0.85,1.2] the PS is naturally divided into three zones which in practice 
are dynamically disconnected due to the presence of invariant KAM spanning curves or their 
remnants. In Fig. we show these zones colored red (zone 1), blue (zone 2) and green (zone 3) 
for k = 0.95. Fet T be the last KAM spanning curve and its point-symmetric partner separating 
zones 1 and 2 (dashed lines in Fig. and the last spanning curve and its point symmetric partner 
separating zones 2 and 3 (solid lines in Fig.[^. 12 is the golden ratio KAM curve. At fc ~ 0.9164.. 
the r and at k = kc the 12 KAM curves are destroyed. For k < 0.91 the term “local” for the 
FFTCFs refers to dynamics within a single zone while the “finite time” is irrelevant since an 
ensemble of chaotic trajectories with initial conditions within a zone remains there for very large 
time scales (~ 10® iterations). Although the chaotic component of the PS becomes connected 
through the breaking of the spanning curves, there is always strong stickiness on the boundaries, 
and, as a result, the mean time needed for a trajectory starting within each one of these zones to 
enter in a neighboring zone, is very large (~ 10"^ — 10® iterations). In fact, even at /c = 1 where 
all spanning curves are broken, a trajectory starting in the red zone needs at least 10 ® iterations to 
enter in the blue zone, while a trajectory starting in the blue zone needs at least ~ 10 ® iterations to 
enter into the green zone. To illustrate this in a more transparent way, we calculate the cumulative 
distribution function Fr(f) = Prob{T < t) of the first passage time T for a chaotic trajectory to 
cross up to time t, either the remnants of T starting from zone 1 , or the remnants of (2 starting 
from zone 2. In practice to find Fr(f) we start with an ensemble of 10® trajectories in a selected 
zone (1 or 2) and we determine the number of trajectories having first passage time T smaller than 
t normalizing by the total number of trajectories in the ensemble. 

In Fig. we show TV(f) for the first passage time from zone 1 to zone 2 (red line) as well 
as the first passage time from zone 2 to zone 3 (blue line). As expected both curves saturate at 
Ft = 1 for f —)■ cx). We clearly observe the lower cut-offs at times ~ 10® and ^ 23 ,mm ~ 10® 

respectively. In addition the corresponding mean times (defined by FT{{t)) = 1/2) are (^ 12 ) = 
2.6 ■ lO"^ and (^ 23 ) = 2.8 ■ 10®. Thus, as dictated by the results shown in Fig. there are well 
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defined time scales which restrict the dynamics to a single zone. This justifies the use of the terms 
“finite time” and “local” in the correlation function of Eq. Q even for chaotic trajectories with k 
larger but close enough to kc, so that the PS stratification into zones, as demonstrated in Fig. 
is still dynamically valid. For a large ensemble of trajectories with initial conditions in the red 
zone, choosing for example N 10^ in Eq. (|3| it is guaranteed that the obtained correlation 
functions characterize mostly the local dynamics within this zone. A similar argument applies 
also for analogous ensembles of trajectories with initial conditions in the other two zones. The 
general conclusion is that FFTCFs are especially useful for the description of correlations at least 
in 2 cases: (i) when the PS of the dynamical system is stratified into zones by invariant spanning 
curves and within each zone local chaos is fully developed (0.85 < k < 0.91 for the SM) and 
(ii) close to the critical point for a dynamical system exhibiting a transition from local to global 
chaos (0.91 < k < 1.2 for the SM). Far beyond the critical point {k 3> kc for the SM), when the 
chaotic sea becomes almost homogeneous and dominating in PS while the remnants of the regular 
dynamics shrink significantly, the FFTCFs converge rapidly, i.e. for relatively small trajectory 
length N, to the usual correlation functions in Eq. Q and become independent of the cell in the 
chaotic domain. 

Since the delay time m should be much smaller than the length of the trajectories involved in 
the calculation of the FFTCFs we have to choose N as large as possible in order to allow for a 
large variation of m. On the other hand, for k > kc, N should be of the order of the minimum 
first passage time in the corresponding zone. For a given N, in order to ensure that the obtained 
numerical results are representative of the limit S'd —?• oo, we use as criterion the convergence of 
2 significant digits when increasing Sd by a factor of 10. In addition the convergence with respect 
to the used length of the trajectories in the ensemble is tested by increasing by a factor of 2. We 
found empirically that taking the values S'rf = 5 ■ 10^ and N = 5 ■ 10^ in the calculation of Q we 
achieve the required convergence. Thus, in the following discussion, if not stated differently, we 
present results using exclusively these parameter values for N and Sd- 


A. Correlation functions in the regular domain 

First we discuss the results for FFTCFs using initial condition cells in the regular domain which 
turn out to possess some universal features. Cells within the regular domain contain quasi-periodic 
and periodic orbits. The corresponding FFTCFs perform oscillations around zero. Increasing Sd 


8 


the amplitude of the oseillations diminishes, tending to zero for large m. This is due to the faet that 
the ensemble averaging is over oseillating forms with slightly shifted frequeneies whieh interfere 
destruetively eaneeling eaeh other. This effeet is demonstrated in Fig. where we plot the LFTACF 
of p in (a) with = 1 and in (b) with Sd = 10®. The LFTCCF in the same eell (not shown here) 
eonverges also to an oseillatory form around zero, deeaying mueh faster (m 30) and having 
initially (m < 30) a very small amplitude (0(10“®)). 

This is a typieal behavior for all eells in the regular domain. Sinee the LFTCFs in the regular 
domain oseillate always around zero we will not eonsider them furthermore in the following. 


B. Auto-correlations in the chaotic domain 

As already mentioned, the LFTCFs may depend on the speeifie loeation of the eell within the 
ehaotie domain. We start our analysis ehoosing as eell the square [0, 0.25] x [0, 0.25] in the 
{6^p) plane whieh is in the immediate neighborhood of the C* order unstable fixed point (0, 0) of 
the SM. We first ealeulate the LFTCFs in this eell and subsequently we explore their dependenee 
on the loeation or the size of the eell. 

In figure 12 ab) we present the LFTACF for k = 0.95 using the phase spaee variables p (a) and 
9 (b) respeetively. Figs. 2^4) show eomponent-wise (p (e), 9 (d)) a typieal trajeetory of the SM 
in zone 1 of the ehaotie domain for this fc-value. 

Let us first diseuss the behavior of the LFTACF^, whieh is mostly small fluetuations around 
zero. This is easily explained by the form of the 6*-time-series. For the eonsidered fc-values 9 per¬ 
forms an oseillatory motion. Consequently also the eorresponding LFTACF^ for a single trajeetory 
possesses an oseillatory form. Averaging over different orbits will induee eaneellations sinee the 
different oseillatory orbits vary in phase, amplitude and frequeney. This is aetually the ease for 
all values of /c, up to fc ~ 1.5 where the orbits beeome strongly ehaotie eharaeterized by irregular 
fluetuations, independently of the eell loeation in PS. Following the same reasoning as it was the 
ease for orbits in the regular domain, we eonelude that the LFTCCFg are fluetuations around zero 
and they also will not be eonsidered furthermore in the following. Thus our subsequent analysis 
foeuses only on the eorrelations of the p-eomponent. 

In eontrast with the LFTACF^, the auto-eorrelation of p, as observed in Fig. |2a) is not trivial 
and the form of the eorresponding trajeetory (Fig.j^c)) is qualitatively different. To explore further 
the auto-eorrelations of p we ealeulated LFTACFp for several values of the parameter k. The results 
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are summarized in Fig. In all ealeulations we used initial eonditions in the eell [0, 0.25] 
[0,0.25]. In eaeh ease the eonvergenee eriteria are well satisfied. Considering the behavior in 
the semi-log seale shown in Fig. [^b), we observe an exponential deeay trend for all values of 
k G [0.6,1.5]. Although the exponential behavior is evident, there are fluetuations of the LFTACFs 
around the exponential envelope. These fluetuations are not of statistieal origin sinee they persist 
by inereasing the number of trajeetories in the ensemble by a faetor 10^. In faet they ean be 
attributed to the stiekiness of the ehaotie trajeetories to the stable manifolds of the set of unstable 
periodie orbits surrounding the stability islands. The LFTACFs deeay to zero signifieantly faster 
for all fc > 1.2 reaehing an almost instantaneous deeay for k > 1.5. The reason for this behavior 
will be diseussed in depth in the next seetion where an orbit analysis is presented. Despite the 
oseillatory fluetuations of the LFTACFs it is useful to perform a fit using an exponential funetion 


ACp{m) = 


(4) 


to determine the eharaeteristie exponent T{k) for the exponential deeay of the auto-eorrelations at 
different fc-values. The resulting T{k) is shown as a funetion of k in Fig.[^ 

For k G [0.6, /Cc] where loeal ehaos is present, the (inverse) exponents reside within a eharae¬ 
teristie plateau, possessing an approximately eonstant value of 650. Above kc where the transi¬ 
tion from loeal to global ehaos takes plaee, the exponents rapidly deerease, approaehing zero for 
k ~ 1.5. In seetion IV we will see that the form of the ehaotie trajeetories ehanges signifieantly 
for k > 1.5 being eharaeterized by irregular fluetuations, while this is not the ease for k < 1.5. It 
will also beeome evident that the transition from loeal to global ehaos induees this ehange in the 
orbits’ strueture, leading to vanishing LFTCFs for k > 1.5. 

Let us now turn to the dependenee of the LFTACFs on the eell’s loeation within the ehaotie 
domain. This is demonstrated in Fig.|7} showing LFTACFs ealeulated for different eell loeations. 
It is evident that the LFTACFs depend on the zone where is loeated. Furthermore, we have 
eheeked the dependenee of these results on the loeation of the eell within eaeh zone. We have 
found that the form of the LFTACF does not ehange as we move within a zone, provided that 
we avoid the mixing of different PS domains (regular, ehaotie). 
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C. Cross-Correlations in the Chaotic Domain 


Proceeding along the same lines as for the case of the auto-correlations, we analyze the LFTC- 
CFs using as the cell [0,0.25] x [0,0.25]. In Fig. [^we present the LFTCCFs with initial 
conditions in this cell for various fc-values. Like in the LFTACFs case, the form of the LFTCCFs 
varies continuously with k and for fc > 1.5 the LFTCCFs approach zero rapidly with increasing 
m (see Fig.j^a)). Despite similarities, the structure of the LFTCCFs is quite different from that 
of the LFTACFs, as it becomes evident if we inspect the LFTCCFs in log-log (Fig.j^b)) and log- 
linear (Fig. [^c)) scales. Notice that the difference of approximately 3 orders of magnitude in the 
initial values of LFTCCFp and LFTACFp relies on the fact that the cross-correlation function of 
each pair of orbits in the corresponding ensemble, is not normalized to one at m = 0 in contrast to 
the auto-correlation functions (see caption of Fig. [^. 

In Figs. [^b,c) it is clearly seen that, neglecting small amplitude fluctuations, the LFTCCFs 
show either a logarithmic or a power-law behavior, thus they possess long-range characteristics. 
In fact for k <kc the LFTCCFs can be better described as a logarithmic function of the delay m, 
while for /Cc < < 1.5 they follow approximately a power-law. Since the long-range character of 

the LFTCCF weakens as k increases beyond k^, it is natural to ask the question whether it is related 
to the local structure of the chaotic domain induced by the last KAM spanning curves. Therefore 
in the following we will try to explain the emergence of long-range correlations by exploring in 
detail the SM chaotic dynamics close to kc. In a recent study llTbll it has been shown that long-range 
cross-correlations can emerge in systems of non-interacting particles, each performing the same 
intermittent dynamics. Thus, one could ask if any kind of intermittent behavior within the chaotic 
domain of the Standard Map for these fc-values is present (at least effectively) and generates these 
long-range cross-correlations. 

Before we address the possibility of intermittency in the SM, a study of the dependence of 
the form of the LFTCCFs on the location of the cell is in order. Concentrating on the case 
k = 0.95 we calculate as a first step the LFTCCFs using cells located in the different zones of 
the chaotic domain as presented by the different colored regions in Fig. The results of the 
corresponding simulations are summarized in Fig.|^ 

In Fig. 1^ it clearly can be seen that the LFTCCFs in zones 2 and 3 have a quite similar form 
while the LFTCCF in zone 1 behaves differently. In the next section, this will be explained by 
exploring the dynamics in the neighborhood of the first order unstable fixed point (0, 0). As in the 
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case of the LFTACFs, the results for LFTCCFs within eaeh zone do not depend on the loeation of 
the eell . However, we observe a very slow eonvergenee when is loeated in the position 
“4” within zone 1 (see Fig. |^. This strange behavior will also be diseussed at the end of the 
next seetion. Notiee that in the ease of zone 2, sinee the ehaotic sea is very thin, one has to use 
partieularly small eells for the ealeulation of the LFTCCFs. Due to the faet that the LFTCCFs in 
zones 2 and 3 are oscillating around zero while the LFTCCF in zone 1 has a long-range profile, in 
the next seetion we will foeus exelusively on the emergenee of eross-eorrelations in this zone. 


IV. INTERMITTENCY AND ORBIT ANALYSIS 

Intermitteney is eharaeterized by long intervals of regular motion (laminar phase) of the trajee- 
tories, interrupted by ehaotie bursts iTTSl . At a first glanee there is no eonneetion of the Eqs. ([T]) 
deseribing the evolution in the SM with the normal form of intermittent dynamies (T] . Therefore it 
is not expeeted to find intermittent eharaeteristies in the first iterate of the SM. On the other hand, 
an effeetive dynamies of the p-variable eould generate sueh a behavior possibly in higher iterates. 
In order to inspeet this hypothesis, a more in depth analysis of the trajeetories in the parametrie re¬ 
gion k E [0.6,1.5], where the loeal to global ehaos transition takes plaee, is required. Furthermore, 
as we have shown in the previous seetion, the long-range slowly deeaying eross-eorrelations whieh 
eould be assoeiated with intermittent dynamies ifT^ oeeur exelusively in SM trajeetories evolving 
within zone 1 of the PS. Thus, in the following we will foeus on dynamieal eharaeteristies of the 
ehaotie trajeetories in this PS region for k < kc- 

In this range of fc-values the p-eomponent of all trajeetories in the ehaotie region of zone 1 
follows the behavior shown in Fig.[T0]^a-e). For k beyond 1.5 this strueture is essentially destroyed 
beeoming gradually a eompletely irregular time-series of random fluetuations. A eareful investi¬ 
gation of the orbit shown in Fig.[T0)^a-e), whieh is representative of all orbits in this ehaotie domain 
for these fc-values, leads to the following eonelusions: 

• The time-series of p has a very intriguing strueture. Considering as “down” the p-values 
below TT and as “up” the p-values above it, we observe (see Figs. [^a,b)) that the p time- 
series eonsists of parts in whieh it is eonseeutively jumping up and down and parts in whieh 
it stays either only up or only down for a large number of iterations. This resembles roughly 
an intermittent time-series where the laminar phase eorresponds to intervals for whieh the 
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p-values are only up or only down while the ehaotie bursts eorrespond to intervals with 
sueeessive up-down (or viee versa) jumps of p. 


Considering the p time-series (Fig. [T^e)) at smaller time seales we observe that it eonsists 
of many sub-intervals eaeh eontaining 4 to 7 points (with p-values all above or all below tt) 
forming an almost triangular strueture. This property is typieal for any p time-series in zone 
1 of the ehaotie domain and it also does not depend on the time interval we zoom into. As 


demonstrated in Fig. 12 it originates from the unstable fixed points of order 4, 5, 6 and 7 of 
the SM. These are densely embedded in zone 1 and their manifolds affeet the trajeetories 
leading to the triangular strueture formation. 


• In faet, the generation of the triangular struetures requires the synergy of the already men¬ 
tioned unstable fixed points and a spanning eurve (or its’ remnants) aeting as a separatrix, 
whieh keeps p eonfined in zone 1. To illustrate this it is useful to eonsider the traee of the 
triangular strueture in PS. Starting from the neighborhood of the PS point at (0, 0) the tra- 
jeetory approaehes the neighborhood of the PS point (27r, 0) in 4 to 7 iterations forming the 
triangular strueture in the p time-series. After reaehing (27r, 0) the trajeetory either returns to 
the neighborhood of (0, 0) in the next iteration, extending the duration of the laminar phase, 
or it jumps to the neighborhood of (27r, 27r) signaling the beginning of a “ehaotie burst”. 
The same proeess ean also take plaee for trajeetories involving the neighborhoods of the 
PS points (27r, 27r) and (0,27r) whieh are the point-symmetrie partners of (0, 0) and (27r, 0). 
During the ehaotie burst interval of the p time-series, the variable 6 performs an oseillation 
from 0 to 27r and baek (full eyele). When the p time-series is in a laminar phase then the 6 
time-series eonsists of half eyeles either from 2tt to 0 (“up” laminar phase of p) or from 0 to 
2ti (“down” laminar phase of p). This behavior is elearly displayed in Fig.fT^f). 

• For k > kc the last spanning eurve is destroyed. Therefore p is not bounded any more and 
the triangular strueture is gradually degraded. However, for /c-values up to 1.5 the rem¬ 
nants of the separatrix still keep for long time intervals the p-trajeetory trapped in zone 1 
and therefore, although slightly deformed, the triangular strueture is preserved. As a eon- 
sequenee the eorresponding LFTCCFs ean exhibit long range eorrelations for delay times 
eomparable with the time seale needed for the p trajeetory to eseape from zone 1. For ex¬ 
ample, at /c = 1.2 although the separatrix is destroyed, orbits with initial eonditions in the 
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neighborhood of (0, 0) and length about N = 10"^ still exhibit the triangular structure with¬ 
out any significant change. However, orbits with N = 10® have escaped through the traps 
created by the residues of the separatrix and contain large parts obeying irregular chaotic 
motion. This behavior is typical for all p-time series when k > kc. 

• The structure of a typical PS trajectory changes completely for k > 1.5, as shown in Fig.fTTj 
Zooming in we observe that the time-series is completely irregular in both 6 and p without 
any structure which could lead to a non-trivial behavior of the LFTCFs. This confirms the 
results presented in the previous section for fc-values in this range. 


The previous discussion suggests to represent the p time-series by defining a suitable symbolic 
dynamics: every triangular unit consisting from 4 to 7 values of p below tt is mapped to “0” and 
every triangular unit consisting from 4 to 7 values of p above tt is mapped to “1”. Obviously the 
opposite choice makes no difference. With this representation of the orbit we can easily define the 
laminar phase as a string of consecutive Is or Os, i.e. 11111111.., and the chaotic bursts as strings 
that have alternating Is and Os, i.e. 10101010... Calculating the LFTCCFs with the symbolic time- 
series yields identical results as those in the previous section, scaled in time by a constant factor 
of ~ 0.2 since every 4-7 points of the original trajectory are mapped to one point in the symbolic 
time-series. 

It is now straightforward to calculate the lengths A of the laminar phases and to obtain the 
laminar length distributions. We used ensembles of 5 ■ 10® trajectories of length N = 1.5 • 10^ 
with initial conditions in the cell [0,0.25] x [0, 0.25], for each fc-value. The results are presented 


in Fig. 13 We checked convergence by increasing the number of trajectories in the ensemble 


by a factor of 10. The result coincides within 3 significant digits with that presented in Fig. 13 


Remarkably enough it turns out that the laminar length distribution as defined above possesses 
some kind of universality. To verify this we recalculated this distribution by changing the size and 
the location of the cell of initial conditions in the chaotic domain of zone 1 and we found, for given 
k, the same result within a two significant digits accuracy. 


The distributions shown in Fig. 13 follow approximately power-laws for k G [0.6, fee]. Es¬ 
pecially at k = kc the laminar length distribution is very close to a power-law. Beyond kc the 
distributions resemble a power-law but this form quickly dissolves, already having deviations at 
k = 1.2. For A: > 1.5 the laminar lengths as well as the associated symbolic dynamics lose their 
meaning, as the orbits no longer posses the aforementioned triangular structure. 
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The fact that for k close to kc the laminar length distributions resemble a power-law is connected 
directly with intermittency. In ifT^ it has been shown that the laminar length distributions for 
intermittent dynamics follow power-laws. There it was argued that the condition for the emergence 
of long range cross-correlations is that the mean laminar length (A) diverges. By applying a power- 
law fit of the form 

PW = ^ ( 5 ) 

to the laminar length distribution dX k = kc, we find a = 1.85. This implies that the average 
laminar length is 

poo poo 

(A) = / P{X)XdX = / A-°-®®dA ^ CX) (6) 

J2 J2 

which is in accordance with the statements in [fT^ . 

Two final comments are in order. As stated in subsection IIIC, in the calculation of the LFTCCF 
within zone 1 we observed a very slow convergence when the cell of initial conditions was 
chosen in the neighborhood of (27r, 0) (cell with index “4” in Fig.|^. Here we will argue that this 
slow convergence is related to the presence of the stable manifold of the first order fixed point in 
this cell. To illustrate this it is useful to derive the corresponding stable and unstable eigenvectors 
in terms of the non-linearity parameter k\ 


-k+s/{‘i: + k)k 

unstable : | ^ 

1 


stable 


—k+^y {4+k)k 


For k ^ 1 these eigenvectors are: 


( 7 ) 


unstable : 



stable 



( 8 ) 


Their direction is shown by the blue lines in Fig.|^ According to this figure the cell with index “1” 
contains the unstable manifold of (0, 0) while the cell “4” contains the stable manifold. Although 
we were not able to determine the exact mechanism which induces the aforementioned slow con¬ 
vergence, we numerically proved that the stable manifold is the origin of this effect. To this end 
we calculated the LFTCCF for the time reversed dynamics (which transform the stable manifold 
to an unstable one) of the SM: 


Hn+l Pn 

Pn+1 =Pn-k sin On+l 


( 9 ) 
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using initial conditions in the cell “4”. The obtained LFTCCF turned out to be identical to that 
using forward in time dynamics and initial conditions in the cell “1”. This clearly verifies our 
previous argumentation. 

The last issue to be discussed concerns the understanding of the qualitatively different behavior 
of the LFTACFs and LFTCCFs for fc-values around kc- For this purpose let us consider more 
carefully the definition of the LFTCFs in Eq. (j^. The sum j^^rCc) in Eq. (j^ implies a summation 
over initial conditions i, j either with i = j (auto-correlations) or with i ^ j (cross-correlations) 
where a;-'^^(0), lie in the PS cell Eocusing on the case i ^ j and assuming that 

ergodicity is valid within a single PS zone, we could consider the chaotic trajectory j as equivalent 
to the trajectory i shifted forward in time by Ur^i iterations, where j is the recurrence time in 
for the trajectory i. Then we could approximately write {n + m) k. {n -f -f m) and the 
sum over initial conditions occurring in the definition of the EETCCEs could be written as: 


=^^P{nr,i) ( 10 ) 


where p{nr^i) is the distribution of the recurrence times in Thus the averaging over different 
trajectory pairs is equivalent to an averaging over single trajectories like in the auto-correlation 
case, performing a random additional time shift following the recurrence time distribution and 
summing over all possible shifts. When the cell is located in the immediate neighborhood of 
the unstable fixed point (0,0) the distribution p{nr^i) can be approximated by the laminar length 
distribution attaining a power-law form. If the cell is small, the distribution p{nr^i) becomes to 

a good approximation independent of i and one can use the notation for Then the EETCE 
in Eq. Q can be written as: 


LCl''\m) = lim tt Pi'^r 

^ Sa^oo Sd ^ ^ ^ 


JV—m—1 


x^f^ {n)xf^ (n + Ur + m) — 


N — m 


n=0 

N—m—1 N—m—1 


(N — m)2 ■4_ 

^ ' n=0 


^ (n) ^ X (n + rir + m) 


n=0 


( 11 ) 


The term in the brackets is the auto-correlation function for the trajectory Xj(n) evaluated at the 
delay time Ur + m. As we have shown in section IIIB the EETACEs in zone 1 follow, for given 
k, an exponential form. The corresponding characteristic exponent T{k) is defined as an ensemble 
property. Calculating the ACEs for individual trajectories forming the ensemble, we find that they 
can be approximated also by exponential functions of the delay m however with an exponent Ti{k) 
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which depends on the initial conditions of the trajectory i. Thus, after performing the summation 
over i we can approximate LFTCCF in zone 1 as: 



( 12 ) 


where p{nr) ~ " with a ~ 1.85 as dictated by the laminar length distribution. The sum over 

can be performed leading to: 



( 13 ) 


with Cn a normalization factor. The sum in Eq. ( [T3] ), as a weighted infinite sum of exponentials, 
leads to the approximate power-law form of the LFTCCF. 

V, CONCLUDING REMARKS 

We have introduced localized finite-time correlation functions as a suitable tool to explore the 
evolution of dynamical systems with mixed phase space. These correlation functions are sensitive 
to local phase space structures, quantifying the impact of stickiness due to phase space traps on 
the emergence of correlations in an ensemble of chaotic trajectories. Using as a prototype model 
the Standard Map we calculated these localized finite-time auto- and cross-correlation functions 
focusing on values of the non-linearity parameter k in the regime of the local to global chaos 
transition. There the phase space is dynamically divided into three zones and the form of the 
aforementioned correlation functions differs in each zone. Specifically, in zone 1 (which includes 
the U* order unstable fixed point) the auto-correlations show an exponentially decaying trend for a 
range of fc-values around the critical kc ~ 0.971635.., while the cross-correlations develop power- 
law tails signaling their long-range character. Around the transition point the dynamics of the 
Standard Map attains intermittent characteristics which can be revealed after a suitable symbolic 
dynamics is introduced. The pathway from this intermittent dynamics to the emergence of long- 
range cross-correlation is similar to the one introduced recently in [fT^ for 1-D dissipative maps of 
Pomeau-Manneville type. Our results demonstrate that intermittency can effectively appear as a 
synergy of complicated phase space networks involving overlaps of unstable and stable manifolds 
of fixed points as well as partial confinement due to remnants of destroyed invariant spanning 
curves. Such a network may in general be established in a low dimensional dynamical system 
close to the critical point associated with the transition from local to global chaos. Furthermore it 
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is confirmed that strong intermittency (even as effeetive dynamies) ean generate long-range eross- 
eorrelations between ehaotie trajeetories as an ensemble property. This seenario ean be easily 
transferred to the emergenee of long-range eross-eorrelations between non-interaeting partieles as 
explained also in [fT^ and may provide a useful manipulation tool for indueing eolleetive behavior 
in non-interaeting systems. 
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FIG. 1. (color online) The mixed PS of the SM for k = 0.95. The colored regions correspond to chaotic 
while the black regions to regular dynamics. The lines indicate the locations of T (dashed) and Q (solid) 
invariant spanning curves. The red colored region defines zone 1 in PS while the blue and the green colored 
regions define zones 2 and 3 respectively (see discussion in section III). 
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FIG. 2. (color online) The cumulative distribution functions for the probability of a chaotic trajectory to 
cross up to time t the second last spanning curve (red line) starting from a PS cell in the first zone (red 
region in Fig. [T]l and the last spanning curve (blue line) starting from a PS cell in the second zone (blue 
region in Fig. [^. Both functions are calculated using k = 1. The dotted line at F'r(f) = 1 is plotted to 
guide the eye. 
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FIG. 3. Numerical demonstration of the cancellation of the correlations in the regular domain. Shown is the 
LFTACF for the p component of orbits in a cell of size 0.5 x 0.5 within the stability island centered around 
(vr, 0) when: (a) a single trajectory, and (b) 10^ trajectories are used in the averaging. The LFTCCFs (not 
shown) have similar form. Analogous results are obtained also for the LFTCFs of the 9 variable. 
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FIG. 4. (color online) Left subfigures depict the average LFTACF for the cell [0,0.25] x [0,0.25] in p ((a)) 
and 6 ((b)) variables for k = 0.95. Right subfigures are typical time-series for p and 6 with initial conditions 
in this cell. Subfigure (c) corresponds to p„ and (d) to 6n versus n. 
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FIG. 5. (color online) (a) LFTACFs for the p variable in the cell [0,0.25] x [0,0.25] for a wide range of k- 
values plotted in linear scale. An ensemble of 5^ = 5 • 10^ orbits each of length N = 5 • 10^ is used, (b) The 
same plot in semi-log scale. The LFTACFp are normalized to one at m = 0 by dividing the auto-correlation 
function of each trajectory in the ensemble with the corresponding standard deviation. 



FIG. 6. (color online) The characteristic exponents T{k) of the LFTACFs (calculated using initial conditions 
in the PS cell [0, 0.25] x [0,0.25]) shown versus k (black circles). The red solid line is an interpolating 
sigmoidal curve to guide the eye. The vertical dotted line indicates the location k = kc- 
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FIG. 7. (color online) The LFTACFs belonging to the different zones of the chaotic domain at fe = 0.95. 
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FIG. 8. (color online) LFTCCFs of p for an ensemble of 5 • 10^ chaotic trajectories with initial conditions 
in = [0, 0.25] X [0, 0.25] (zone 1) obtained using k = 0.8 (black line), k = kc (red line), k = 1.2 (blue 
line) and k = 1.5 (green line), (a) Linear-linear, (b) linear-log and (c) log-log scale. The cross-correlation 
function of each pair of trajectories in the ensemble is normalized by dividing with the root of the product 
of the standard deviations of those trajectories. 
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FIG. 9. (color online) The LFTCCFs calculated for k = 0.95 using PS cells located in different zones of the 
chaotic domain. For zone 1 we use = [0,0.25] x [0, 0.25], for zone = [0.02,0.17] X [1.66,1.8] 
and for zone 3 = [4.82, 5.03] x [3.57, 3.68]. Each cell is numbered with a violet-colored number. The 

blue lines above cells 1 and 4 correspond to the fixed point’s (0,0) eigenvectors (see end of section IV). 
The blue line above cell 1 represents the unstable eigenvector whereas the blue line above cell 4 represents 
the stable eigenvector. 
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FIG. 10. (color online) Time-series of the variables p (a,b,c) and 9 (d,e,f) of the Standard Map, for a typical 
orbit within the chaotic domain, for k = 0.95 and initial conditions (0.1,0.2). Shown are (a,d) the first 10“^ 
iterations, (b,e) a zoom-in in the region [0,10^] and (c,f) a zoom-in in the region [400, 500] of this orbit. 
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FIG. 11. (color online) A typical trajectory for k = 1.8 and initial conditions (0.2, 0.1). In the zoom-in it is 
clearly seen that the triangular substructure in the time-series of p is not sustained and a laminar phase (as 
defined in section IV) does not exist. 
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FIG. 12. (color online) Fixed points of the Standard Map for A; = 1.0 found by the method introduced in 

m. 



FIG. 13. (color online) Laminar length distributions p{X) for various values of k, shown in log-log scale. We 
used 5-10® trajectories with initial conditions in the cell [0,0.25] x [0,0.25] and verified the independence 
of the cell location in zone 1. 
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